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ABSTRACT 

In this paper we examine the problems of phasing using light curves and offer 
an alternate technique using the changes in acceleration to establish the zero point. 
We give astrophysical justification as to why this technique is useful and apply the 
technique to a selection of Type II Cepheids. We then examine some limitations of 
the technique which qualify its use. 
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1 INTRODUCTION 

When studying intrinsically pulsating stars a key challenge 
is to find common physical processes at a given pulsational 
phase in stars of different spectral characteristics and pe- 
riods. In order to do this, a common analysis technique is 
required which can find links between apparently disparate 
groups of objects. This is especially true when trying to 
understand what the stars are physically doing by linking 
spectroscopic features to the photometry. Uncertainties ex- 
ist when using the photometry to establish the zero point, 
as it cannot be guaranteed that the same spectroscopic fea- 
tures are occurring at the same photometric phase between 
different stars. This leads to ambiguities in comparisons be- 
tween stars. An alternative approach is required that allows 
a clearer comparison and a better understanding of the phys- 
ical processes involved. 

In this paper we examine the problems of phasing us- 
ing light curves. We offer an alternative technique using the 
minimum of the acceleration curve as the zero point. We 
give astrophysical justification as to why this technique is 
useful and apply the technique to a selection of stars. Some 
of the limitations of the technique are then examined. 



2 TRADITIONAL PHASING 

In order to be able to understand the regular pulsations of 
stars one needs to characterize the basic parameters of what 
changes over a pulsation cycle. One method is to plot the 
data (such as photometry or radial velocity) as a function of 
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pulsational phase. This is shown in Figure^for a selection of 
Type II Cepheids, with periods ranging from 1.6 d (SW Tau) 
to 90 d (U Mon). In this, using a fixed zero point and a 
known period, data over a wide time interval can be com- 
pared. Traditionally, the zero point of this phasing has been 
based on the time of maximum (e.g. lMoffett fc Barnesll984l : 
ISchmidt et al.l2003allb|) . or minimurn lpollard et alll99d) of 
the star's light curve. This is due to the historical precedent 
that stellar variability is generally first discovered through 
overall light variations. In the more regular Cepheids such 
as SW Tau and k Pav (top, Figure 0, the determination of 
the position of an extremum is a relatively straight-forward 
process. However the more complicated RV Tauri stars (i.e. 
U Mon, lowest curve, Figure^, are less obvious. There is an 
additional problem that the phases of maximum and min- 
imum light do not necessarily occur at the same time in 
different photometric filters. As such, while the extrema in 
the light curves represent a relatively easily observed quan- 
tity, they do not represent an optimal means by which to 
compare the astrophysical behaviour of the stars. Two main 
problems are encountered in phasing data based on maxi- 
mum or minimum light. Firstly, a photometric filter region 
must be chosen, and secondly, the light extrema must be 
well defined. 

The first is trivial only if all the photometry used is 
from a single band-pass. As soon as other band-passes are 
observed, the extrema do not exactly coincide because they 
are sampling different regions of the stellar atmosphere. This 
is clearly seen in Figures EK an d Et an d Table Q where 
maximum light shifts to a later phase at longer wavelengths 
(from B to 7). Astrophysically, this is due to changes in tem- 
perature and hence flux distribution over the pulsation cycle. 
As the stellar photosphere expands and cools, the peak flux 
shifts to longer wavelengths, causing maximum light in the 
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Figure 1. Photometry (left) and velocity curves (right) of selected Type II Ccpheids phased using a zero point of minimum V magnitude. 



redder band-passes to occur slightly later in the pulsation 
cycle. An extreme case is W Vir (Figure and Table 
in which maximum light shifts forward in phase by 0.2 of a 
pulsation cycle between V and R as the primary maximum 
becomes less prominent in the red. This maximum is most 
likely associated with the shock wave propagat ing through 
the star during this phase (lAbt fc Hardid ll960). 

The second problem, an ill-defined maximum or min- 
imum, will depend critically on the data density obtained 
and the repeatability of the light curve. As with any fitting 
to observational data, a large scatter in the data points can 
introduce large errors in the extrema fitting. The broadness 
of the extrema can also add to the difficulty in determining 
the zero point. This is seen in SW Tau (Figure HJl where the 
minima are quite broad in all the band-passes observed. As 
seen from Figure Q this represents one of the more regular 
stars of the sample. The scatter for U Mon is even greater 
and is affected by long term variations in the light curve and 
irregularities in the depths of the minima. 

Both of these problems lead to slight shifts in the zero 
points between the stars, such that the same phase does not 
represent the same physical state for each star. It should 



also be noted that different stars spend a different fraction 
of their pulsation cycle brightening from minimum to max- 
imum light. Hence, when basing a zero point on minimum 
light, the phases of maximum light compared are offset and 
vice versa. This can been seen in Figures ISK and0Ji and 
Table where the duration of the rise from minimum to 
maximum light ranges from 0.2 to 0.5 of a pulsation cycle. 
This leads to difficulty in comparing stars with a range of 
light curves at the same phases. 

An alternative approach which phases the astronomi- 
cal data based on the same physical state for each star is 
desirable. 



3 ALTERNATE PHASING TECHNIQUES 

The pulsation cycle of a star can be observed physically in 
a number of ways through changes in brightness (photom- 
etry), radial velocity (spectroscopy) and angular diameter 
(interferometry) . Having discussed the problems of using the 
light curves, we turn to the other techniques. Both give an 
indication of radii changes, though the angular diameter far 
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Table 1. Phase shifts" (columns 2 to 9) of photometric extrema relative to minimum V light, which by definition is phase=0.0. 
Columns 10 and 11 are respectively the phase shift between photometric and dynamic zero points (A) and its uncertainty (HWHM). 
See text (end of Section for detailed explanation of these terms. 
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"The uncertainty in these phase shifts varies from star to star and between passbands. The range is from ±0.07 (for V in W Vir) 
to ±0.16 (for / in re Pav) when measured at 0.1 mag above the minimum (or below the maximum) of any light curve. 
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Figure 2. SW Tau (a) B VRI photometry, (b) spline-fitted mean 
Fe I radial velocities, and (c) acceleratio n curves. S W Tau is 
phased on a 1.58356-day period from|McSavcncy ( 2003]) and using 
the zero point from the acceleration curve. 
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Figure 3. n Pav (a) BVRI photometry, (b) spline-fitted mean 
Fe I radial velocities, and (c ) acceleration curv es, re Pav is phased 
on a 9.0714-day period from McSavcne^ 120031) and using the zero 
point from the acceleration curve. 



more directly than the radial velocities. Maximum, mean 
or minimum radius would be useful as a calibration point. 
To date Type II Cepheids have not been observed interfer- 
ometrically, so the radial velocities remain the most easily 
observed parameter. Ideally, one would use radial velocity 



curves to obtain radii changes and hence well denned zero 
points for phasing the data. However, this is not a straight- 
forward process. 

Derivation of radii changes from velocity curves re- 
quires integration of the pulsational velocity curve. This 
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Figure 4. W Vir (a) BVRI photometry, (b) spline-fitted mean 
Fe I radial velocities, (c) acceleration curves, and (d) He* and He I 
equivalent widths . W Vir is phased on a 17.2768-day period from 
iMcSav cncv ( 2003) an d using the zero point from the acceleration 
curve. 



differs from the measured radial velocity curve due to limb- 
darkening effects across the disk and the fact that the star is 
a spheroidal object changing in size and not a point source 
moving towards and away from us. This can be corrected 
for by the use of a geometric projection and limb dark- 
ening coefficient p, but the exact value of t h e coe fficient 
varies in the litera t ure (1 .31, Iffindslev fe Bel3 <198fjft . 1.41, 
lAlbrow fe Cottrelll (119941) ). As well as varying from star to 
star, this "constant" also varies as a function of pulsational 
phase, method of velocity measurement and lin e formation 
depth JSabbev et alJll995l : iNardetto et alll2004l) . These val- 
ues and statements are also based on discussions of Classical 
Cepheids and their applicability to the Type II Cepheids is 
not proven. In addition, the radius derived from a particu- 
lar line or species will trace the radius of that particular line 
formation region. Certain infra-red regions will probe deeper 
into the stellar photosphere and produce different radii to 
those determined from iron lines in the visual wavelength 
regions. 

Minimu m radius of the star has been used previously 
for phasing (iMarengo et al.ll2002ft by setting the zero phase 
point to be the time at which the pulsational velocity moved 
from a positive to a negative value. This does, however, re- 
quire the transformation of the radial velocity into a pulsa- 
tional velocity, with t he uncertainties ment ioned earlier. In 
a more recent paper, IMarengo et al] |2004) advocated dy- 
namical phasing based on the radial velocity curve. The rea- 
soning given was this phasing was based on hydrodynamical 
quantities and could thus be more readily compared with 
the time-dependent hydrodynamical models. These authors 
also discussed the phase lag observed using both the light 
maximum and radius minimum zero points. 

Given the provisos of working with the pulsational ve- 
locities and radii, it was decided to examine the acceleration 
curve produced by taking the derivative of the radial veloc- 
ity curve. While not the true acceleration curve, as it is not 
the derivative of the pulsational velocity, it is a clear indi- 
cation of changes in the bulk motion of the line formation 
region. The minimum of this acceleration was then used, as 
for these stars it marked the transition between the bulk 
of the material in the line formation region falling inwards 
and moving outwards. By using the same species from the 
same wavelength regions, a more robust comparison of ma- 
terial under the same temperature and pressure conditions 
is obtained for the different stars. This also minimized the 
number of assumptions required in analyzing the data and 
is model independent. 

In applying this technique a spline curve is fitted to the 
phased radial velocity curve of a particular line or species 
(see Figures^, |3J) and^fc). Care is taken to fit over several 
repeated cycles to avoid end effects. This curve is then differ- 
entiated to find the acceleration. The minimum of this curve 
defines the phasing zero point, namely minimum accelera- 
tion. The data are then rephased based upon this. Figures 
[5J:, and 0J; show that the acceleration minimum is far 
more clearly defined than the photometric minima shown in 
Figures 03., EK and Hi- 

Quantitative estimates of the uncertainty of the derived 
dynamical zero point phase have been determined from the 
acceleration data, based upon the half width at half the 
minimum of the acceleration dip (see HWHM in Table 1). 
We note that it is not just dependent upon the data density 
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around the minimum acceleration point, but also depends 
upon the steepness of the radial velocity curve at that point. 



4 APPLICATION 

To test this technique of dynamical phasing, it has been 
applied to a selection of Type II Cepheids. The observa- 
tions are taken from McSaveney (2003). The BVRI pho- 
tometry is from observations made by A.C. Gilmore and 
P.M. Kilmartin, using automated, single-channel photome- 
ters mounted on the Optical Craftsmen 0.61-m telescope at 
Mount John University Observatory (MJUO). The veloci- 
ties are measured from spectra obtained by JAM, using a 
high-resolution spectrograph mounted on the McLellan 1.0- 
m telescope at MJUO. Spectroscopic data were reduced us- 
ing FIGARO processes which were supplemented by MAT- 
LAB scripts. 

The velocity curves shown here were obtained by aver- 
aging a selection of Fe I line velocities. These Fe I lines, with 
a range of excitation potentials and line strengths, represent 
a region of the star's photosphere in which line formation 
is well understood and can be considered representative of 
the star. We are aware o f line level effects betwee n various 
species (see for example IWallerstein et all (I1992T) ). but in 
this paper we are considering only those layers representing 
these Fe I lines. 

The individual velocities were measured by taking a 
mean of the line bisector velocities at depths 0.7, 0.8 and 0.9 
for each line. Th i s follo wed the line bisector method as deve l- 
oped bv lAlbrowl Jl994f) and used in lWallerstein et al.Nl992l) . 
This was used in preference to methods using line core ve- 
locities, to smooth out any line asymmetries. These line 
asym metries were outside the scope of this work jMcSavenevI 
l2003h . 

The three stars shown in this paper (SW Tau - Figure 
[5] ft Pav - Figure [3] and W Vir - Figure constitute a 
representative sample of Type II Cepheids of periods ranging 
from 1.58 to 17.28 days. The Fe I lines from which these 
velocities were measured show a range of behaviour over a 
pulsational cycle, from little change in line width in SW Tau 
(the shortest period star), to mild broadening at phase 0.0 
in k Pav, to extensive broadening and line-splitting at phase 
0.0 in W Vir. As can be seen from Figures EU and 01 the 
dynamical phasing technique can be used on all of these 
stars, despite the differences in velocity curve produced by 
the different line behaviours. 

The acceleration curves shown for each star (Figures^, 
and^JJi) all have distinct minima that have been rephased 
to 0.0. Even in the case of SW Tau, where a second minimum 
is found around phase 0.4-0.5 due to a bump in the velocity 
curve, this is easily distinguished from the primary velocity 
field reversal. 

The result of this phasing provides a much easier com- 
parison between these stars, as it is clear when the motions 
of the velocity fields are reversing. In contrast, by phasing 
using V minimum (Figure 0, it occurs between phase 0.85 
and 0.95 in SW Tau, at ~0.8 for k Pav and ~0.7 in W Vir. 
These are shown quantitatively in Table 1, where the phase 
shifts between V minimum and dynamical zero points (A) 
are given for the five Type II Cepheids in our sample. These 
As show that there is more than 0.2 in phase difference 



between these values that is not related in any way to a 
physical stellar parameter. Also the dynamical zero points 
have a much better defined position (HWHM in Table 1 is 
between 0.01 and 0.10) when compared to the uncertainties 
in the photometric zero points, which range from 0.07 to at 
least 0.16 and vary considerably from passband to passband 
for a given star. This effect is also clearly visible in Figures 
[5] [3] and 2] where the width (in phase) of the acceleration 
minimum for a given object is in all cases much narrower 
than any of its respective photometric passband minima. 

Astrophysically, this shows that the dynamical zero 
point provides a clear common phase compared to that indi- 
cated by any particular photometric bandpass. The dynam- 
ical zero point indicates when the direction of motion of the 
line formation region shifts from falling inwards to moving 
outwards. In these stars, this is when shock-wave features 
are detected as a shock propagates through the line form- 
ing regions. Hence spectral features such as line broadening, 
splitting and emission all occur around this phase. This is 
seen in the Ha and He I equivalent widths (see Figure |ljl) 
where the strongest emission occurs at around phase 0.0. 
This is also observed in a wider range of Type II Cepheids 
(McSaveney, 2003). 



5 DISCUSSION 

There are several assumptions which will have an affect on 
the results obtained when using this technique. These in- 
clude what line species are used, the method of measuring 
the radial velocities and how the spline curve is fitted. 

An important point to note with this technique is that 
different line species are formed over slightly different re- 
gions of the stellar photosphere, and will therefore give 
slightly different zero points. While this may be construed 
as a problem as large as that of the shifts in zero point from 
the photometry, this is not actually the case. The spectral 
line formation regions can be far more specifically defined 
than the far broader photometric regions, particularly for 
the dominant species (usually Fe I), from which the veloc- 
ities are determined. When establishing the velocities care 
should be taken to use lines of similar species, excitation 
potential and \og(gf) such that there is a high degree of con- 
fidence that the lines are formed in the same region of the 
stellar photosphere. Averaging Ha velocities with Fe I ve- 
locities would be counterproductive given the large phase 
lags and amplitude differences obse rved between these lines 
in Type I and Type II Cepheids j jWah^sl^nej^T] Il992l : 
IVinko et alJll998t |Pettersonll2002t iMcSavenevlbOfm 

The radial velocity measuring technique will also affect 
the velocity determined. Techniques fitting to the sides of 
the lines, or taking the line bisector value from midway up 
or higher (closer to the star's continuum) in the line, will be 
more easily distorted by line broadening or splitting. This 
will result in different amplitudes. Techniques fitting to the 
line bisector from deeper in the line core (as used in this 
work) give a better indication of the maximum and mini- 
mum velocities achieved by material in the line forming re- 
gions. The crucial point is to be consistent in the manner 
of application such that it is clear that the stars (as shown 
by the appropriate spectral lines) show the same measured 
zero-point behaviour. 
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As with any data set, good radial velocity phase cover- 
age is essential, particularly at the critical phases between 
minimum and maximum light when the minimum radius oc- 
curs. To avoid distortion by end effects in the spline fitting, 
repetition of the full velocity curve is required. This tech- 
nique is quite data intensive. However, for the comparison 
of stars with a range of periods all with large datasets, it is 
invaluable. 



6 CONCLUSION 

Using the acceleration changes as a zero point in dynamical 
phasing offers a substantial improvement over previous phas- 
ing techniques in understanding the pulsations in Type II 
Cepheids. It shifts the phasing emphasis from the more am- 
biguous light extrema to the astrophysically more specific, 
yet still observable, changes observed through radial velocity 
curves. This allows a better understanding of the physics of 
the stellar pulsations and clearer comparison between stars, 
as they can be collectively phased to a more clearly defined 
astrophysical state. 

This technique has been applied to a selection of Type II 
Cepheids and found to work extremely well for a variety of 
velocity curve shapes and periods. From a preliminary sur- 
vey of SW Tau, k Pav and W Vir, it is clear that the shock- 
associated features of line broadening and line splitting, and 
Ha and He I emission, are associated with the bulk motion 
reversal observed and can be more clearly compared by us- 
ing this particular zero point. More detailed study of the 
shock behaviour for a wider range of Type II Cepheids us- 
ing this technique can be found in McSaveney et al. (2005, 
in preparation). 
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